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Abstract. We report on key studies on the dynamics of black holes (BHs) in gas-rich galaxy 
mergers that underscore the vital role played by gas dissipation in promoting BH inspiral down 
to the smallest scales ever probed with use of high-resolution numerical simulations. In major 
mergers, the BHs sink rapidly under the action of gas-dynamical friction while orbiting inside 
the massive nuclear disc resulting from the merger. The BHs then bind and form a Keplerian 
binary on a scale of < 5 pc. In minor mergers, BH pairing proceeds down to the minimum scale 
explored of 10-100 pc only when the gas fraction in the less massive galaxy is comparatively 
large to avoid its tidal and/or ram pressure disruption and the wandering of the light BH in 
the periphery of the main halo. Binary BHs enter the gravitational wave dominated inspiral 
only when their relative distance is typically of ~ 10~'^pc. If the gas preserves the degree of 
dissipation expected in a star-burst environment, binary decay continues down to O.lpc, the 
smallest length-scale ever attained. Stalling versus hardening below < O.lpc is still matter of 
deep investigations, and there is no unique answer depending on the yet unexplored dynamics 
of gas in the vicinity of the binary. 



1. Introduction 

Dormant black holes (BHs) with masses in excess of ^ lO^M© are found to be ubiquitous in 
bright spheroids, today [H [2]. This local population comprises the dead remnants of a bright 
past, when the same BHs were powering the most luminous quasars. The recent discovery of 
scale relations between the BH mass and the properties of the stellar bulge [3], likely set via 
AGN feedback, has prompted the study of BH growth and evolution in the general framework 
of galaxy evolution. According to the current ACDM paradigm for structure formation, galaxies 
interact and merge as their dark matter halos assemble in a hierarchical fashion [3], and BHs 
incorporated through mergers into larger and more massive systems evolve concordantly: major 
central gas-inflows are triggered during the violence of a merger that feed the BH and power 
AGN activity [5j. In this context, close BH pairs form as inescapable outcome of galaxy evolution 
[6]. In our local universe, NGC 6240 and Mrk 463 provide compelling evidence of ongoing gas- 
rich mergers where two active nuclei are present, still at large separations of ~ kpc. Whether 
these BHs will spiral inward, form a BH binary and coalesce under the emission of gravitational 
waves is a matter of our concern. The Laser Interferometer Space Antenna (LISA) is expected 
to record these extraordinary events out to redshift z ~ 20 providing not only a firm test 



of General Relativity, but also a view, albeit indirect, of galaxy clustering together with an 
extremely accurate measure of the BH mass and spin [71 [U [9] . 

Galaxy mergers cover cosmological volumes (of hundred kpc aside) whereas BH coalescences 
probe volumes from a few parsecs (when they bind in nuclear discs) down to an astronomical unit 
and less. Thus, following a merger, how can BHs reach the gravitational wave inspiral regime? 
Our aim is at studying the BH dynamics in gas-rich environments, and in particular, the transit 
from the state P of pairing when each BH moves individually inside the time-varying potential 
of the colliding galaxies, to state B when the two BHs dynamically couple their motion to form 
a binary. After all transient inflows have subsided and the new galaxy has formed, the BH 
binary, surrounded by a massive circum-nuclear disc, enters phase H where it hardens to smaller 
separations under the action of gas-dynamical and gravitational torques, ideally down to ~ 10~^ 
pc (for typical BH mass of W^Mq) where the gravitational waves domain G starts. There is a 
number of key questions to address: (i) How does transition from state P ^ B depend on the 
gas thermodynamics and level of dissipation? (ii) In the grand nuclear disc inside the remnant 
galaxy, how do orbits evolve? (iii) Do the BHs reach the gravitational wave driven domain? (iv) 
During the hardening through phase B — > H, do the BHs collect substantial amounts of gas to 
form cold individual discs? 

2. Pairing of Massive Black Holes in gas-rich mergers 

There are two types of mergers: major mergers between galaxies of comparable mass (1:1 mass 
ratio), and minor mergers between galaxies with smaller mass ratios (1:10 typically). 

Major Mergers has been studied with N-Body/SPH simulations with unprecedented force 
resolution (down to ~ 1 pc using splitting techniques and GASOLINE as integrator; see [LQl) to 
describe the collision of two galaxies similar to the Milky Way. Each galaxy comprises a central 
BH of 2.6 X 10^ Mq, a stellar bulge, a disc of stars and gas (with mass fraction of 10% relative 
to the total disc mass), and an extended spherical dark matter halo (of 10^^ Mq) with NFW 
density profile (see [TU] for details). 

The galaxies first experience two close fly-bys: in this early phase, the cuspy potentials of 
both galaxies are deep enough to allow for the survival of their baryonic cores that sink under 
the action of dynamical friction against the dark matter background, dragging together the two 
BHs. Strong spiral patterns appear in both the stellar and gaseous discs, and as the merger 
continues, non-axisymmetric torques redistribute angular momentum: as much as 60% of the 
gas originally present in each disc of the parent galaxies is funneled inside the inner few hundred 
parsecs of the individual cores. This is illustrated in the upper right panel of Figure 1, where 
the enlarged color coded density map of the gas is shown, after 5.1 Gyr from the onset of the 
collision. Each BH is surrounded by a rotating stellar and gaseous disc of mass ~ 4 x 10^ Mq 
and size of a few hundred parsecs. The two discs and BHs are just 6 kpc far apart, and at the 
same time a star-burst of ~ SOM© yr~^ has invested the central region of the ongoing merger. 

At this stage, the simulation is stopped and restarts with increased resolution (of ~ 2 pc). In 
order to simulate the environment of a star burst where cool gas coexists with the warm phase 
heated by stellar feedback, the pressure is set equal to P = (7 — l)pu with 7 = 7/5 (according 
to fits by [llj). The internal energy per particle u evolves with time as a result of PdV work 
and shock heating modeled via the standard Monaghan artificial viscosity term. 

With time, the two baryonic discs get closer and closer and eventually merge in a single 
massive self- gravitating, rotationally supported nuclear disc, now weighing 3 x 10^ Mq. This is 
illustrated again in Figure 1 (mid and bottom right panels). The gaseous disc, dominant in 
mass, is surrounded by a background of dark matter and stars distributed in a spheroid. 

The BHs have been dragged together toward the dynamical center of the merging galaxies, 
and move inside the grand disc spiraling inward under the action of gas-dynamical friction. In 
less than a million years after the merger, they eventually bind gravitationally to each other. 




Figure 1. The different stages of the merger between two identical disc galaxies seen face-on. 
The color-coded density maps of the gas component are shown using a logarithmic scale, with 
brighter colors for higher densities. The four panels to the left show the large-scale evolution 
at different times (obtained with a force resolution of 100 pc). The boxes are 120 kpc on a 
side (top) and 60 kpc on a side (bottom) and the density ranges between 10~^ atoms cm~'^ and 
10^ atoms cm~'^. During the interaction, tidal forces tear the galactic discs apart, generating 
spectacular tidal tails and plumes. The upper panel to the right shows a zoom in view of the 
two discs before they merge into a single rotating nuclear gaseous disc embedded in a series 
of large-scale ring-like structures (middle panel). The boxes are now 8 kpc on a side and the 
density ranges between 10~^ atoms cm~^ and 10^ atoms cm~'^. The two bottom panels, with a 
gray color scale, show the detail of the inner 160 pc of the middle panel (here the force resolution 
is 2 pc); the nuclear disc is shown edge-on (left) and face-on (right), and the two BHs are also 
shown in the face-on image. 

as the mass of the gas enclosed within their separation is less than the mass of the binary. It 
is the gas that controls the orbital decay, not the stars. The transition between state P to B 
is now completed as illustrated in Figure 2. Dynamical friction against the stellar background 
would bring the two BHs this close only on a longer timescale, ~ 10* yr [TU]. This short sinking 
timescale comes from the combination of the fact that gas densities are much higher than stellar 
densities in the center, and that in the mildly supersonic regime the drag against a gaseous 
background is stronger than that in a stellar background with the same density [12] • It is worth 
noticing that the transition P ^ B is sensitive to the gas thermodynamics: BH coupling is 
delayed if gas were to follow thermal evolution with a 7 = 5/3 [lOj . 

Does BH pairing proceed similarly, in minor mergers predicted to be common events in the 
high redshift universe [8j and of primary importance for LISA [9]? To answer this question 
we extended our numerical investigation to 1:4 mergers at z = 0, and 1:10 mergers at 2; = 3, 
assuming a roughly constant Mbh — -^buige relation in between these cosmic epochs, and initial 
galaxy models replica of a Milky Way suitably rescaled in mass and size (see for details [T3] . 
The masses of the two BHs in the z = 3 runs are thus 6 x 10^ and 6 x 10^ Mq, and their expected 




Figure 2. Orbital separation of the two BHs as a function of time during tlie last stage of the 
galaxy merger. The orbit of the pair is eccentric until the end of the simulation. The two peaks 
at scales of tens of parsecs at around t = 5.1213 Gyr mark the end of the phase during which the 
two holes are still embedded in two distinct gaseous cores. Until this point the orbit is the result 
of the relative motion of the cores combined with the relative motion of each BH relative to the 
surrounding core, explaining the presence of more than one orbital frequency. The inset shows 
the details of the last part of the orbital evolution, which takes place in the nuclear disc arising 
from the merger of the two cores. The binary stops shrinking when the separation approaches 
the force resolution limit (2 pc). 



inspiral and coalescence signal falls nicely in the LISA sensitivity window [H]. 

It is found that minor mergers differ profoundly from major mergers as early noticed by [14j . 
The encounter is closer to an accretion process whereby the less massive galaxy is dramatically 
damaged during its sinking into the primary. In our recent study [T3] , paring is found to be very 
sensitive to the details of the physical processes involved. In all cases with no-gas (i.e., in "dry" 
runs) the formation of a close BH pair is aborted: tidal shocks progressively lower the density in 
the satellite until it dissolves, leaving a wandering black hole in the remnant. This is illustrated 
in Figure 3 (thin lines) where the BH relative distance remains as large as 1-10 kpc. Only with 
the inclusion of a cold gaseous disc component, and star formation the outcome of the merger 
changes significantly. 

Figure 3 depicts the the stellar and gaseous components of the satellite to show their profound 
structural damage (the primary is not shown). For mass ratios 1:4 at z = 0, bar instabilities 
excited at pericentric passages funnel gas (present in a fraction /gas of the total disc mass) to 
the center of the satellite, steepening its potential well and allowing its survival against tidal 
disruption down to the center of the merger remnant. As shown in Figure 3 (thick lines), the 
BHs pair down to ~ 100 pc scales (the force resolution limit), creating conditions favorable to 
the formation of a BH binary. The smaller satellites (with 1:10 mass ratio at z = 3) are more 
strongly affected by both internal star formation and the gas-dynamical interaction between 
their interstellar medium and that of the primary galaxy. Torques in the early stages of the 
merger are not acting to concentrate gas to the center, due to the absence of a stellar bar and 
the stabilizing effect of turbulence. As a result, ram pressure strips all of the ISM of the satellite. 




Figure 3. BH separation as a function of time in four of our simulations. Upper row: BH 
distance in 1:4 mergers (for galaxy models at z = 0); the thin and thick lines refer to simulations 
with no-gas and with gas (/gas = 0.1%), respectively. Lower row: BH distance for the 1:10 
mergers (for galaxy models at z = 3); the thin and thick lines refer to simulations with no-gas 
and with gas (/gas = 0.3%), respectively. The insets show the color-coded density maps of stars 
(left) and gas (right), 4 kpc on a side. The large dot on the BH curve indicates the time at 
which the two snapshots are recorded. Colors code the range 10~^ — 1 Mq pc~'^ for stars, and 
10~^ — 0~^ Mq pc~^ for the gas. These snapshots are representative of the average behavior of 
the satellites during the first two orbits. Note the formation of a strong bar for the 1:4 minor 
merger, which is absent for the 1:10 case, and the truncation of the gaseous disc in the 1:10 
satellite caused by ram pressure stripping. 

Only the gas-rich satellites (those with /gas = 0.3) undergo a central burst of star formation 
during the first orbits which increases their central stellar density allowing for their survival. In 
these models, pairing of the two BHs via dynamical friction occurs down to ~ 100 pc, a few 
Gyrs after the disruption of the satellite (see Figure 3). 

3. Black hole binaries in massive nuclear discs 

As shown in Section 2, massive cir cum- nuclear discs form in the aftermath of a major gas-rich 
merger. It is in these discs that the BHs complete their transition from P— > B, and continue 
to spiral inward under the action of gas-dynamical torques from B — > H [161 [T71 118j . A still 
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Figure 4. BH separation as a function of time. One BH is set at the center of a rotationally 
supported disc while the secondary BH moves initially on an eccentric (e = 0.7) either co- or 
counter-rotating orbit. BH masses are of 10^ Mq,. The gaseous disc of mass Mjisc = 10^ is 
described by a Mestel profile. Solid (blue) line refers to the co-rotating case, while (red) dashed 
line refers to the counter-rotating case (Dotti et al. in preparation). 

open issue is whether the BHs will reach the domain of gravitational waves inspiral within a 
Hubble time: for a 10^ Mq BH binary on a circular orbit, the transition H — > G occurs when 
the separation is around 10~'^ pc. Can material and gravitational torques be effective in driving 
the BHs down to this tiny scale? 

Here we describe our attempts to explore the transition from P ^ B ^ H, for BHs orbiting 
inside a circum-nuclear disc using GADGET as N-Body/SPH code with a force resolution of 
only ~ 0.1 pc, and no splitting during the entire course of the evolution [15]. In our selected 
model, a BH (called primary) is set at the center of a massive differentially rotating gaseous disc 
in equilibrium with a stellar bulge (see [T71[18] for the details). A second BH (called secondary) 
with similar/equal mass is delivered at a large distance (~ 50 pc) from the center, along a 
coplanar eccentric (e = 0.7) orbit that can either be co- or counter-rotating relative to the 
background disc. The 10^ SPH particles, making the disc, evolve as in [lOj: accordingly, the 
gas thermodynamics is described by the index 7 = 7/5 that accounts for the presence of net 
cooling in a star-forming region. Shocks here are less important as the equilibrium disc is only 
mildly perturbed by the BHs. Figure 4 shows the BH separation as a function of time for the 
co-rotating and counter-rotating cases. No stalling is observed in both cases, as the relative BH 
distance decays rapidly (in < 20 Myr) down to the force resolution length-scale. In the final 
stages, the more rapid orbital decay is due to the torque exerted by the ellipsoidal deformation 
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Figure 5. z-component of the orbital angular momentum of the secondary BH normalized to its 
initial value Lq for the counter-rotating case, described also in Figure 4, to show the occurrence 
of the angular momentum flip induced by dynamical friction [15] . 

that forms in the gas when the heads of the density wakes overlap and whose axis is misaligned 
relative to the BH binary axis |16j . 

In the counter-rotating case, the angular momentum of the secondary BH (initially negative) 
grows very efficiently during the first Myr, when the BH is passing through the central, high 
density region of the disc. Angular momentum continues to grow monotonically for the next 
3 — 4 Myrs, then becomes positive, i.e. the BH starts to move on a co-rotating orbit with respect 
to the disc. This is illustrated in Figure 5. It is dynamical friction that causes this orbital 
"angular momentum flip" . In both cases the eccentricity of the orbit decreases to very small 
values (also in the counter-rotating case, after the orbit becomes co-rotating) due to the different 
response of the fluid to the gravitational pull of the BH at the different orbital phases. When at 
pericenter the BH moves faster than the gas and it is decelerated by the density wake excited 
behind its trail; when at apocenter the BH is moving more slowly than the gas and the wake is 
trailing in front causing a tangential acceleration. The composite effect is a decrease of e. 

In a suite of runs, the BHs have been modeled as "sink particles", i.e. they are allowed 
to accrete gas particles during their dynamical evolution [15]. We introduced an "on flight" 
algorithm for accretion and determined the amount of gas that binds to the BHs. It is only 
when the BH binary circularizes that gas is accreted to such an extent that both BHs are 
surrounded by their own accretion disc, and these discs are expected to play a role in guiding 
the subsequent hardening phase down to the gravitational domain. 

4. Open issues 

Numerical simulations, carried on with unprecedented accuracy, have revealed that the transit 
of dual BHs from P — > B ^ H, and finally from H^ G is a sensitive function of the merger 
type and of amount of cold gas present in the interacting galaxies. While the transition from P 



(pairing) — > B (binary formation) appears to be likely in gas-rich major mergers as well as in gas- 
rich minor mergers at high redshift, hardening down to the gravitational wave domain remains 
still uncertain on scales below ~ 0.1 pc and not fully explored. It has been suggested that a 
circum-binary viscous disc inevitably forms around the BH binary on sub-parsec scales that 
absorbs the the angular momentum of the binary [19j . This circum-binary disc would represent 
the last cold environment for BH hardening from H — > G. Braking of the BH rapid motion 
requires energy loss and angular momentum transport trough a mechanism that is reminiscent 
of planet migration in proto-stellar discs [19; ,20j : while tidal torques from the BH binary carry 
away orbital angular momentum, viscous torques inside the disc sustain the radial motion of 
the gas toward the BHs, maintaining the binary in near contact with the disc. Equilibrium 
between these two torques would cause the slow drift of the BHs toward smaller and smaller 
separations, until gravitational waves guide the final inspiral. No calculation has reproduced yet 
the formation of a circum-binary disc from the earlier phase, in a self-consistent manner, nor 
it is clear how fast will be the inspiral, and how large the growth of the eccentricity [21]. The 
BH binary likely enters phase G with a residual eccentricity still imprinted in the gravitational 
wave signal, despite the circularizing action of the gravitational wave back reaction. LISA is 
expected to be lunched by 2020. By that time, our theoretical understanding of binary hardening 
in a gas-rich environment will hopefully improve thanks to the progress expected in numerical 
simulations, and in our ability to model fragmentation, star formation and feedback, inside 
galactic nuclei. 
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